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1. Introduction. 



Understanding of low energy phenomena in QCD, such as confinement and chiral 
symmetry breaking or, in more general terms, the strong coupling problem and a ground 
state structure in an asymptotically free non-abelian gauge theory is, without doubt, one 
of the main (if not the main) problems in modern quantum field theory. In spite of years 
of attempts to answer this question we are still far from complete satisfaction, although 
a lot of interesting and promising ideas were suggested during the first 20 years of QCD 

i- 

Considerable progress has been made in this direction during last years using numer- 
ical approach of the lattice gauge theory 0. The lattice gauge theory calculations are, 
however, still incomplete. Apart from that, they sometimes leave behind an unpleasant 
aftertaste (although this is of course, a very subjective matter) that one obtains numer- 
ical results without gaining real understanding of the underlying physics. To our minds, 
understanding of these issues in the framework of an analytical approach would be in- 
valuable. An analytic method that is capable of solving the low energy sector of QCD 
starting from first principles, would also, presumably, teach us a lot about other strongly 
interacting theories such as technicolor. 

Unfortunately, the arsenal of nonperturbative methods to tackle strongly interacting 
continuum theories is very limited, to say the least. Methods that perform very well 
in simple quantum mechanical problems, are much more difficult to use in quantum 
field theory (QFT). This is true, for example, for a variational approach. In quantum 
mechanics it is usually enough to know a few simple qualitative features in order to set 
up a variational ansatz which gives pretty accurate results, not only for the energy of a 
ground state, but also for various other vacuum expectation values (VEV). In QFT one 
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immediatelly is faced with several difficult problems when trying to apply this method, 
as discussed insightfully by Feynman @. 

First, there is the problem of calculability. That is, even if one had a very good guess 
at the form of the vacuum wave functional (or, for that matter, even knew its exact 
form) one would still have to evaluate expectation values of various operators in this 
state. In a field theory in d spatial dimensions this involves performing a d- dimensional 
path integral, a problem, in itself very complicated and, in general not manageable. This 
problem is especially severe in nonabelian gauge theories, where gauge invariance poses 
strong restrictions on admissible trial wave functionals (WF). In this case it becomes very 
difficult to find a set of WF's which are both, gauge invariant and amenable to analytic 
calculation. 

Another serious problem is the problem of "ultraviolet modes" . This means the follow- 
ing. In a variational calculation of the kind we have in mind, one is mostly interested in 
the information about the low momentum modes. However, the VEV of the energy (and 
all other intensive quantities) is dominated entirely by contributions of high momentum 
fluctuations, for a simple reason, that there are infinitely more UV modes than modes 
with low momentum. Therefore, even if one has a very good idea, how the WF at low 
momenta should look like, if the UV part of the trial state is even slightly incorrect the 
minimization of energy may lead to absurd results. Due to the interaction between the 
high and low momentum modes, the IR variational parameters will in general be driven to 
values which minimize the interaction energy, and have nothing to do with the dynamics 
of the low momentum modes themselves. 

Even though over the years many attempts at variational calculations in QCD have 
been made these two problems invariably made their presence felt, and at this point 
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one really can not point to any succesful variational calculation in a nonabelian gauge 
theory. Our feeling is, however, that these obstacles are not necessarily insurmountable, 
and that this direction is still far from being exhausted and deserves further development. 

In this paper we present a variational calculation of the vacuum WF in a pure SU(N) 
Yang Mills theory, which, at least partially, is free from the problems mentioned earlier. 
We use wave functionals that are explicitly gauge invariant. The correct UV behaviour 
is built into our ansatz. In the case at hand this can be done due to the assymptotic 
freedom of the models considered. We are able to calculate VEVs of local operators in 
our trial states in a reasonable approximation, combining the renormalization group and 
the mean field techniques. 

Our main result is that the energy is minimized at the value of the variational pa- 
rameter away from the perturbative vacuum state. This leads to a dynamical generation 
of scale in the vacuum WF. The value of the vacuum condensate - < F°;„F":„ > in this 
state turns out to be equal to O.OOSGef^ for Aqcd = 150Mev. Even though this result 
should be taken only as an order of magnitude estimate (due to approximations made), 
it is pleasing to see a number so close to the phenomenologically known 0.012Gev'^ |^ 
emerge from this simple calculation. 

The paper is organized as follows. In Section 2 we set up the variational calculation and 
discuss in some detail our variational ansatz. In Section 3 we discuss the approximation 
scheme for calculation of VEVs in the trial WF. In Section 4 the minimization of energy 
and calculation of < > are performed. Some elements of the calculation of the Wilson 
loop and an area law is discussed in Section 5. Finally, Section 6 contains dicussion of 
our results and outlines directions for future work. 
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2. The variational trial state and the gauge invariance. 



As mentioned in the previous section, an immediate question one is faced with, when 
picking a possible variational state is calculability. One should be able to calculate aver- 
ages of local operators in this state 

<0>= J D(/)^*[(/)]0^[(/)] (2.1) 

A calculation of this kind obviously, is tantamount to evaluation of a Euclidean path 
integral with the square of the WF playing the role of the partition function. One should 
therefore be able to solve exactly a c?- dimensional field theory with the action 

^[0] = -logvl/*[0]vl/[0] (2.2) 

Since in dimension d > 1 the only theories one can solve exactly are free field theories, 
the requirement of calculability almost unavoidably restricts the possible form of the WF 
to a Gaussian (or as it is sometimes called squeezed) state: 

m = exp |-i I Sxd'y [<P{x) - C(x)] G-\x, y) [<P{y) - ((y)]} (2.3) 

with C(x) and G{x,y) being c-number functions. The requirement of translational invari- 
ance usually gives further restrictions: ({x) = const, G{x,y) = G{x — y). 

The restriction to Gaussian WF is of course a severe one. However, one can still hope 
that in some cases the simple Gaussian form can capture the most important nonpertur- 
bative characteristics of the true vacuum. Indeed, the Gaussian variational approximation 
has been used succesfuUy in self interacting scalar theories, where it is known to be exact 
in the limit of large number of fields. Perhaps the most celebrated use of these trial states 
is the BCS calculation of the superconducting ground state 0, where for most of the 
interesting quantities its accuracy is of order 10-20%. 
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The reason the approximation works well in these theories, is that in both cases a 
single condensate dominates the nonperturbative physics, and the Gaussian ansatz is 
wide enough to accomodate this most important condensate. From this point of view, it 
would seem then, that it is perfectly reasonable to try a similar variational ansatz in the 
Yang - Mills theory. After all, it is strongly suggested by the QCD sum rules that the 
pure glue sector is dominated by one nonperturbative condensate < >. We also know, 
that the VEV of the field strength itself < F > vanishes, since it is not a gauge invariant 
operator. A state of the form ( p.3| ) with C = would indeed give zero classical fields, but 
nonzero quadratic condensates. 

There is, however, one obvious difficulty with this idea. It is very easy to see, that 
in a nonabelian theory it is impossible to write down a Gaussian WF which satisfies the 
constraint of gauge invariance. The SU(N) gauge theory is described by a Hamiltonian 



H = d^x 



(2.4) 



where 



Efix) = I ^ 



5AHx) 



Btix) = ^e.,,{d,At{x)-d,A^{x)+gr"^A%x)Al{x)} (2.5) 
and all physical states must satisfy the constraint of gauge invariance 

G^{x)'$[A] = [diE^ix) - gr^^A^i{x)E^{x)] ^[A] = (2.6) 
Under a gauge transformation U (generated by G'"(x)) the vector potential transforms as 

A^{x) ^ Af"(x) = S^\x)A^{x) + \nx) (2.7) 



where 



S''\x) = hr (t'^U^t^'u) ; A^(x) = -tr (r'^U^d^u) (2.8) 
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and r'* are traceless Hermitian N by N matrices satisfying tr[T°'T^) = 26°'^. A gaussian 
wave functional 



1 



nA^] = exp |-- J d'xd'y [A1{x) - Q{x)] [A]{y) - C^{y)\^ (2.9) 

transforms under the gauge transformation as 

^[A^] ^ vI/[(A^)«] (2.10) 

In the abehan case it is enough to take diG~j^ = to satisfy the constraint of gauge 
invariance, In the nonabelian case, however, due to the homogeneous piece in the gauge 
transformation ( |2.7| ), no gauge invariant Gaussian WF exist. 

One possible strategy is to disregard this fact and hope that one does not loose much 
by minimizing the energy in the whole Hilbert space, which also includes unphysical states. 
This is, however, very risky. The sticking point is that the hamiltonian of the theory is 
unique only on physical states. One can add to equation (|2.4|) an arbitrary operator 
multiplied by one of the generators of the gauge group without changing the energies of 
the physical states, but reshuffling the rest of the spectrum beyond recognition. In this 
way the gap between the physical vacuum and some of the unphysical states can be made 
very small. In fact the large Hilbert space can even contain states which have energies 
lower than the physical vacuum. Since we are working with a particular Hamiltonian, it is 
not clear a priori that this is not the case. Therefore minimizing the energy on the whole 
space may lead to huge admixtures of unphysical states in the "best variational state", 
making the results of such a procedure meaningless. Of course, one could be lucky and 



with the particular choice of the Hamiltonian (|2.4| ) all unphysical states may have large 



energies, but there is no way to know it without a separate investigation of this question. 

We, at any rate, will restrict our attention to gauge invariant states only. It is clear 
then, that the Gaussian ansatz must be modified. Several modifications were considered 



in pevious work. One obvious possibiity is to restrict classical fields to zero and insert 
adjoint Wilson lines in the exponential 0, so that 

- Q{^)] iG-')t^{x, y) [A]{y) - C}(y)] - Bt{x)G~\x - y)B]{y)W^\C) (2.11) 

where W{C) = Pexp {ig Jq dliF°-A^), and F°- are the generators of SU(N) in the adjoint 
representation. This form, however, makes it practically impossible to perform explicit 
calculations, except in the week coupling limit. Another proposed modification is to 
multiply the Gaussian by a finite order polynomial in the fields. In that way gauge 
invariance can be maintain to a finite order in the coupling constant 0. Then, however, 
it is again not quite clear to which extent the calculation is nonperturbative. 

Instead, we will take a straightforward approach, and simply project the Gaussian WF 
onto gauge invariant sector. In this paper we also restrict ourselves to the case of zero 
classical fields (C = 0). Our variational ansatz is therefore 

^A^] = J DU{x)exp |-1 J d'xd'y Af^{x)Gr/-'{x - y) Af{y)^ (2.12) 

with A^"- defined in (|2.7| ) and the integration is performed over the space of special unitary 
matrices with the SU (N) group invariant measure. 

Before attempting a calculation with this expression, we will impose several restrictions 
on the form of G, which will lead to considerable simplifications. First, we will only 
consider matrices G of the form 

Gl!;{x-y)=5'''5,,Gix-y) (2.13) 

This form is certainly the right one in the perturbative regime. In the leading order in 
perturbation theory, the nonabelian character of the gauge group is not important, and 
the integration in equation ( |2.12| ) is basically over the U{1)^^~^ group. The S""^ structure 
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is then obvious - there is a complete democracy between different components of the 
vector potential. The 5ij structure arises in the following way. If not for the integration 
over the group, would be precisely the (equal time) propagator of the electric field. 
However due to the integration over the group, the actual propagator is the transverse 
part of G^^. It is easy to check that the longitudinal part diG^j^ drops out of all physical 
quantities. At the perturbative level, therefore, one can take Gij ~ 6ij without any loss 
of generality. We will adopt this form of the matrix G also in our variational calculation. 

We can use additional perturbative information to restrict the form of G even further. 
The theory of interest is asymptotically free. This means that the short distance asymp- 
totics of correlation functions must be the same as in the perturbation theory. Since 
G~^ in perturbation theory is directly related to correlation functions of gauge invariant 
quantities (e.g. E'^), we conclude 

G-\x)^\, x^O (2.14) 
x^ 

Finally, we expect the theory nonperturbatively to have a gap. In other words, the 
correlation functions should decay to zero at some distance scale 

G(x)~0, x>^ (2.15) 

We will build this into our variational ansatz in the simplest possible way. We will take M 
to be our only variational parameter. This can be done by choosing for G{x) a particular 
form that has the UV and IR asymptotics given by (|2.14| ) and ( p.l5|) , like, for example a 



massive scalar propagator with mass M. We find another parametrization slightly more 
convenient. The form that will be used throughout this calculation has the following 
Fourier transform: 
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We have checked, that using a massive propagator instead, practically does not change 
the results. Equation ( p.l2|) together with equations (|2.13|) and ( p.l6|) define our varia- 
tional ansatz. We now have to calculate the energy expectation value in these states and 
minimize it with respect to the only variational parameter left - the scale M. Note, that 
the perturbative vacuum is included in this set of states, and corresponds to M = 0. A 
nonzero result for M would therefore mean a nonperturbative dynamical scale generation 
in the Yang - Mills vacuum. In the next section we will explain the approximation scheme 
we use to calculate expectation values in the trial state. 



3. The effective sigma model and the renormaUzation group. 



The question now is, how do we calculate expectation values in the state of the form 



<0> = ^ J DUDU' <0>A (3.1) 

^0>A = I DAe'^^ / dxdyAY-ix)G-\x-y)AY-iy) q J dx' dy' A^'\x')G-\x' -y')AV'\y') 



where Z is the norm of the trial state. Two simplifications are immediatelly obvious. 
First, since we will only be considering gauge invariant operators O, one of the group 
integrations is redundant. Performing the change of variables A A^ (and remembering 
that both inegration measures DU and DA are group invariant), we obtain (omitting the 
volume of SU{N) factor / dU) 

<0> = ^ f DU <0 >A (3.2) 



Z 

^0>A = I DAe'^ / dxdyAf-(x)G-Hx-y)AY-iy)Q^-^, J dx' dy' A){x')G'Hx' -y')A](y') 
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Also, since the gauge transform of a vector potential is a linear function of A ( |2.7|) , for 
fixed U{x) tliis is a Gaussian integral, and can therefore be performed explicitly for any 
reasonable operator O. We are left then only with a path integral over one group variable 
U{x). But this is a tough one! 

Let us consider first the normalization factor Z. After integrating over the vector 
potential we obtain: 



Z = J DUexp{-T[U]} (3.3) 

with an action 

T[U] = ^TrlnM + ^X[G + SGS^]-^X (3.4) 

where multiplication is understood as the matrix multiplication with indices: colour a, 
space i and position (the values of space coordinates) x, i.e. 

{AB)r,{x,z) = J d'yA^^{x,y)B%{y,z), XOX = J d'x<fyXt{x)0'i{x-y)X]{y) (3.5) 



The trace Tr is understood as a trace over all three types of indices. In equation (|3l^) we 
have defined 

St^{x, y) = S'"'{x)6i,6{x - y), M^^{x, y) = [S^^'{x)S'\y) + 6^']G-\x - y)5,, (3.6) 

where S^^x) = \tr {t'^UW'U) and A^(a;) = Hr (r^U^diU) were defined in (U) and tr 
is a trace over colour indices only. Using the completeness condition for SU (N) 

Tt-rt, = 2 (6a6,, - (3.7) 



ij'kl ^ y^it^jk 

one can see that S""^ is an orthogonal matrix 
where we used that tr [Ut^W) = trr" = 0. 
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We have written action (|3.4| ) in a form which suggests a convenient way of thinking 
about the problem. The path integral ( p.3|) defines a partition function of a nonlinear 
sigma model with the target space SU{N)/Zjq in three dimesnional Euclidean space. The 
fact that the target space is SU{N)/Zn rather than SU{N), follows from the observation 
that the action ( p.4|) is invariant under local transformations belonging to the center of 
SU (N). This can be trivially traced back to invariance of A" under gauge transformations 
that belong to the center of the gauge group. 

We note, that the quantity U{x) has a well defined gauge invariant meaning, and it 
is not itself a matrix of a gauge transformation. A contribution of a given U{x) to the 
partition function eq. p.3| and to other expectation values corresponds to the contirbution 
to the same quantity from the off diagonal matrix element between the initial Gaussian 
and the Gaussian gauge rotated by U{x). Therefore, if matrices U{x) which are far from 
unity give significant contribution to the partition function, it means that the off diagonal 
contribution is large, and therefore that the simpleminded non gauge invariant Gaussian 
approximation (which neglects the off diagonal elements) misses important physics. 

The action of this sigma model is rather complicated. It is a nonlocal and a nonpoly- 
nomial functional of U{x). There are however two observations, that will help us devise 
an approximation scheme to deal with the problem. First, remembering that the bare 
coupling constant of the Yang Mills theory is small, let us see how does it enter the sigma 
model action. It is easy to see, that the only place it enters is the second term in the 
action ( |3.4] ), because A^(x) has an explicit factor 1/g. Moreover, it enters in the same 
way as a coupling constant in a standard sigma model action. We can therefore easily set 
up a perturbation theory in our sigma model. With the standard parametrization 

U{x) = exp{i^(/.V"} (3.9) 
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one gets Af (x) = —di(f)°'{x) + 0{g), S"''^{x) = 5"* + 0{g) and the leading order term in 
the action becomes: 

/ d^xd^ydi^''{x)G-\x - y)di^''{y) (3.10) 



16 

This is just a free theory, except that the propagator is nonstandard, and at large momenta 
its Fourier transform behaves like 

D{k) G{k)^^ -L (3.11) 

Nevertheless, the perturbation theory is straightforward. Indeed, it is easy to see, that in 
this sigma model perturbation theory the coupling constant renormalizes logarithmically. 
The first order diagram that contributes to the coupling constant renormalization is the 
tadpole. In a sigma model with a standard kinetic term this diagram diverges linearly 
as / d^k/k'^, a sign of perturbative nonrenormalizability. In our model, though, due to a 



nonstandard form of the kinetic term ( 3.11 ), the diagram diverges only logarithmically 



as / d^k/k^. The form of the (3 function therefore is very similar to the (3 function in the 
ordinary QCD perturbation theory. In this paper we assume, that to one loop, the two (3 
functions indeed coinside. The explicit calculation in the framework of the sigma model 
will be presented elsewhere |jlO|. The perturbation theory, therefore becomes worse and 
worse as we go to lower momenta, and at some point becomes inapplicable. 

Now, however, let's look at the other side of the coin. Let us see how does the 
action look like for the matrices U{x), which are slowly varying in space. Due to the 
short rangedness of G{x), clearly for U{x) which contain only momenta lower than the 
variational scale M the action is local. In fact, with our ansatz ( |2.16| ) it becomes the 
standard action 

TlIU] = ^tr J d^x diU\x)d,U{x) + ... (3.12) 
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where we omit the higher order in g terms. We have used the completeness condition 
(131 ) and the fact that tr{WdiU) = to rewrite 

X^{x)X'^{x) = -{l/g^)tr {r^U^diU) tr (r^U^diU) = -{2/g^)tr {u^diU U^d.u) (3.13) 

In this low-momentum approximation we also neglected the space dependence of SfHx) 



gets SGS^ G. 

Strictly speaking, due to the local symmetry of the original theory (^.41), the action 
for the low momentum modes is slightly different. The derivatives should be understood 
as Zjv covariant derivatives. The most convenient way to write this action, would be to 
understand U{x) as belonging to U{N) rather than SU{N) and introduce a U{1) gauge 
field by: 



This defines a sigma model on the target space U {N)/U{1), which is isomorfic to SU {N)/Z]^ 
This action does not look too bad. Even though it still can not be solved exactly, it is 
amenable to analysis by standard methods, such as the mean field approximation, which 
in 3 dimensions and for large number of fields should give reliable results. 

The suggestion therefore, is the following. Let us integrate perturbatively the high 
momentum modes of the field U{x). This is the renormalization group (RG) transforma- 
tion. We would like to integrate out all modes with momenta k"^ > M^. This procedure 
will necessarilly generate a local effective action for the low momentum modes. At the 
same time, because of the (presumable) equivalence of the RG flows in QCD and our 
effective sigma-model, the effective coupling constant will be the the running QCD cou- 
pling constant aQcoiM) at scale M. This part of the theory can then be solved in the 
mean field approximation. Clearly, in order for the perturbative RG transformation to 



in the term SGS'^ 



in (|3.4| ), then using the fact that S is an orthogonal matrix ( |3.8|) one 




(3.14) 
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be justified, the QCD running coupling constant at the scale M must be small enough. 
Our procedure will then make sense, provided the energy will be minimized at the value 
of the variational parameter, for which 

«qcd(M) < 1 (3.15) 

We will check whether this consistency condition is satisfied at the end of the calculation. 
In the next section we will calculate the expectation value of the Hamiltonian in the lowest 
order of this approximation scheme, and perform the minimization with respect to M. 

Before doing that, we would like to make one side remark. It is amusing to see how 
the present framework can accomodate instanton effects. Recall, that in a path inte- 
gral formalism instantons describe the tunneling transition between some initial state 
and a new state ^[A\ where a field Ai is obtained from a field Ai by a large 
gauge transformation which is described by a nontrivial element of the homotopy group 
Yiz{SU {N) / Z n) = Z. The target space of the effective sigma model SU{N)/Zn, has the 
right topology: Il3{SU (N) / Z n) = Z. The model therefore must have classical "hedge- 
hog" solutions analogous to Skyrmions [|1T|. In fact in the perturbative regime they should 



be easy to find. At week coupling the action reduces to (up to a numerical coefficient) 
/ cPxd^ytr W {x)diU{x) (3.^^)4 {y)diU {y) , and equation of motion for U (x) becomes rel- 
atively simple. Note also, that this action has a dilatational invariance x Ax, Skyrmion 
solution must approach 1 asymptotically at large distances. These functions Uci{x) then 
correspond to contributions of the off diagonal matrix elements between the initial Gaus- 
sian, and the same Gaussian gauge transformed by a large gauge transformation, which 
is precisely the meaning of one instanton contribution to the path integral. 

Note that the dilatational invariance is broken in our ansatz for slowly varying modes, 
by apearance of the scale M. Indeed, the only Skyrmion solutions in the low momentum 

15 



effective action (|3.14|) , are pointlike, due to Derrick's collapse. This means physically, 
that the scale M sets the nonperturbative infrared limit on the instanton size. Our 
variational vacuum therefore, is free from the infrared problem associated with the large 
size instantons. 

The variational ansatz which has been considered corresponds to a zero value of the 
QCD 9 - parameter, since we have integrated over the entire gauge group without any 
extra phases. As is well known, the general 6'-vacuum is defined as 

|^>=^e™V> (3-16) 

n 

where n labels the topological sectors in the configuraion space (space of all potentials 
y4^(a;)). Generalization of our trial wave functions to nonzero 9 is trivial - all we need to 



do is to insert in equation ( p.l2|) an extra phase factor in the integrand 



exp J dxeijktr [(f/t9,f/)(f/^9,f/)(f/^4[/)] | (3.17) 

The integrand here is a properly normalized topological charge, and it takes integer values 
for topologicaly nontrivial configurations U{x), i.e. this factor reproduces the exp(m^) 
term in ( ^.161 ). This phase factor can be obtained also if one remembers that usually the 
6'-dependence of the wave functional is given by the ex'p[i9Scs{A)], where Scsi^) is a 
Chern-Simons term, which under the gauge transformation U transforms as 

ScsiA"") = Scs{A) + dxe.,,ktr [{U^ d,U){U^ d,U){U^ duU)] (3.18) 



thus integrating over U leads precisely to the phase factor (|3.17| ). The state thus con- 
structed, is an eigenstate of an operator of the large gauge transformation with eigenvalue 
e*^. This will result in addition of the same topological term to the effective action ( p.4|) . 
It is amusing that for 9 = n, the "Skyrmions" in the effective theory will be "fermions". 
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In the rest of this paper, we shall ignore instanton contributions, but it will be interesting 
to come back to this question later. 



4. Solving the variational equation. 

We will now calculate the expectation value of the energy. 

+ \ge.,keu^r'' j d'xd.AlA'^A^^ + ^e.jkeiimr'T'' J d'xA'^AlAfA^^ (4.1) 

We first perform the Gaussian integrals over the vector potential at fixed U{x). Let 
us consider, for example, the calculation of the chromoelectric energy: 

/ £x < Ef >A= i dPx < - , s , ^ >A = 

J ' J 5Al{x)5Al{x) 

TrG-^ - J SxSyd^zG'\x - y)G"\x ~ z) < A'l{y)A1{z) >a (4.2) 

Using ( |3.4| ) it is easy to calculate the average over A. Defining for convenience 

<(x) = / dSd'z\\{y)G-\y - z)S'\z){M-'r{z, x) (4.3) 



so that gaussian integration over A is / DAexp [—{\/2){A + a)M.[A + a)] one gets 

j d^x < E^ 2 >A= 3{N^ -1) J d^xG~\x,x) - J d^x{G-^M-^G-^)'l^{x,x) (4.4) 

- 1 d^xd^ya1{x)G-\x-y)a1{y) 

where G~^{x—y) = J d^zG~^{x—z)G^^{z—y) andA^"^ is defined as / d^yM~^{x,y)M{y, z) 
6^{x — z). Let us note that G^^ has dimension [x\~^ and M.^^ has dimensiion For 
chromomagnetic field the calculations are straightforward and one gets 

< {e,,,d,Alf >A= {e,,ud,alf + e,,,eamd^d^{M-')Tm{^.y)l=y (4-5) 
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<djAlA'iAl>A = djata'ia^^ + djat{M-')t{x,x) (4.6) 
+ al'aj(A^-^)-(a;,|/)U=, + a;;,aj(A1-i)^^(a;,|/)U=, 

ei^keumr'^'r''' < A'^AlAfAl 2r''^r''^ay,afal (4.7) 
+8r''T'^a'laf{M-'y'{x, x) + 12^^" P^^M'^f^x, x){M-^Y\x, x) 

Here we have used the obvious notation M-l^ = M."'''dij The next step is to decompose 
the matrix field U (x) into low and high momentum modes. In general this is a nontrivial 
problem. However, since we are only going to integrate over the high momenta in the 
lowest order in perturbation theory, for the purposes of our calculation we can write 

U{x) = Ul{x)Uh{x) (4.8) 

where Ul contains only modes with momenta k'^ < M^, and Uh has the form Uh — 
1 + igT"'(f)'}j and (f)H contains only momenta k'^ > M^. This decomposition is convenient, 
since it preserves the group structure. Also, since the measure DU is group invariant, we 
can write it as DUlDUh- With this decomposition we have: 

Xt{x) = S'^^{x)XUx)-rXU^) (4.9) 

Further simplifications arise, since we only have to keep the leading piece in cjf^. We can 
therefore write in our approximation: 

S''\x) = Sf{x) 

M''\x,y) = 25''^G-\x-y) 

X^x) = XUx) + XtH{x) 

= lxt^{x) + ^X'lH{x)Sl^{x) (4.10) 
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We are now in the position to rewrite different pieces in the VEV of energy in this 
approximation: 



I d'xd'yXUx)G-'ix - y)XUy) -I J d'xd'yXU^)G-\x - y)A^^(y) (4.11) 

The cross term vanishes, since to this order, as we shall see , there is a decoupling between 
the high and the low momentum modes in the action, and therefore the product factorizes, 
and < X°'fj >= 0. Our ansatz for G"^ ( p. 16 ) allows us to simplify this expression further. 
Remember that Xl{x) contains only momenta below M. Then it is immediate to see that, 

' d^xd'yXtLG'\x-y)X1M = M^ j dxXIMXt^x) (4.12) 

We can then rewrite equation ([4.11| ) as 

/ d'x < Et ' >A= J G-\x,x) - 

/\/f2 /■ 1 r 

— j d'xXt^ix)Xt^ix) - - y d'xd'yXlHix)G~\x - y)XtM (4-13) 

The contribution of the magnetic term to the energy is very simple. All cross terms 
between the low and high momentum modes drop out. Some vanish for the same reason 



as the cross terms in equation ([4.11|) , and others because they are explicitly multiplied 
by a power of the coupling constant. Since our approximation is the lowest order in g, 
except for the nonanalytic contributions that come from the low mode effective action, 
those terms do not contribute. In fact, the entire low momentum mode contribution 
drops out of this term. The reason is that the only terms which could give a leading order 
contribution, is 

/ ie.jkd,Xt^)' (4.14) 



It can be rewritten as 

' lijL 
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ift^Lf + Oig') (4.15) 



Where f°ji^ is the "magnetic field" corresponding to "vector potential" Af^,- However, 
has the form of a pure gauge vector potential. Therefore /""jj^ = 0, and the contribution of 
this term is higher order in g"^. We have checked, that including this term, indeed changes 
the energy density in the best variational state by a small amount (0(10%)), but has no 
effect at all on the best value of the variational parameter M. The entire magnetic field 

contribution to the energy is then: 

1 1 AT^ — 1 

2 < >A= -^{e,,,d,\lHf + 9r9fG(x - y)l=y (4.16) 

The last step is to perform an averaging over the ?7-field. For convenience, we rewrite 
here the complete expression for the energy density (here V = J cPx is a space volume): 

-G-^(x,x) + (iV2 - l)afafG(x -y)|, 



<2ff > _ 3(iV^ 



- 17 y ^'"^ < A^Jx)A^^(x) >u (4.17) 
where the averaging over the U - field should be performed with the sigma model action 

(|3.4| ). In our approximation this action has a simple form. Using equation (|4.10| ) we 
obtain 

1 r Mr 

T = -j dxdyKH{x)G-\x - y)KH{y) + rfxA^^(x)A^^(x) (4.18) 

The low momentum mode part is precisely equal to in equation (|3.14|) . The only 



difference is, that the coupling constant that appears in this action should be understood 
as the running coupling constant at the scale M. This, obviously is the only 0(0) effect 
of the high momentum modes on the low momentum effective action. 

Tl = -^^^tr [ d^x{di - iAi)U^{x){di + iAi)U{x) (4.19) 



2g\M) 

We are now in a position to evaluate the VEV of energy. The contribution of the 
high momentum modes is immediatelly calculable. Using the parametrization Uh{x) = 
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1 — ^g(f)"'T"', we find that 0" are free fields with the propagator 



\y) >=25-'[d^d\G-\x-y)]-' 



(4.20) 



Also to this order Xini^) = di(j)°'{x) and therefore ^ijkdjXkH = 0- Using (|4.2CI| ) one can see 
that 

1 r N"^ — I d^k 

- J d'xd'y < X1h{x)G-\x - y)XUy) >u= J^,G-\k) (4.21) 

where A is the ultraviolet cutoff, and the the contribution of the high momentum modes 



to the energy (first two lines in equation ( |4.17|) ) is: 

2^0 .,.o .A d^k 



V 



{N' - 1) 



- 1 r fM 



27r2 



(27r)3 
k'^dk 



G-\k) + PG{k) 

-M + ^ +2 / edk\ 

2 M\ Jm J 



+ U -^G-\k) 



2 7o (27r)3 



lOvr^ 



M^ + ... (4.22) 



Terms denoted by ... in eq. [4.22| depend on A, but are independent of the variational scale 
M. 

We now have to evaluate the contribution of the low momentum modes. It is clear 

from the form of the action ( 4.19| ), that this contribution as a function of M will not 

be featureless. The most convenient way to think about it, is from the point of view 

of classical statistical mechanics. Comparing equations ( [4.17|) and ([4.19 ), we see that 

we have to evaluate the internal energy of the sigma model (with the UV cutoff M) at 

the temperature proportional to the running coupling constant g'^{M). For large MQ, 

the coupling constant is small, which corresponds to the low temperature regime of the 

sigma model. In this regime the global SU (N) ® SU (N) symmetry group of the model 

is spontaneously broken. Lowering M, we raise g'^{M), and therefore the temperature. 

At some critical value gc, the model will undergo a phase transition into the unbroken 
* Large M, of course means large relative to Aqqd- 
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(disordered) phase. Clearly, in the vicinity of the phase transition all thermodynamical 
quantities will vary rapidly, and therefore this is a potentially interesting region of coupling 
constants. 

Before analysing the phase transition region let us calculate E{M) for large M. In this 
regime the low momentum theory is weekly coupled. The calculation is straightforward, 
and to lowest order in gives: 

1 AT^ — 1 

-M' < XUx)XtL{x) >= -J^M' (4.23) 



Putting this together with the high momentum contribution, we find 

-M^ M»Aqcd (4.24) 



V 1207r2 

This indeed is the expected result. The energy density monotonically increases as M^, 
with the slope which is given by the standard perturbative expression. Note, however that 
the slope is very small, and the contribution of the low momentum modes to the energy 
is negative. Therefore, if the internal energy of the sigma model grows significantly in the 
phase transition region, the sign of E{M) could be reversed]] and the energy will then be 
minimized for M in this region. 

To see, whether this indeed happens, we will now study the low momentum sigma 
model in the mean field approximation. We rewrite the partition function by introducing 
a (hermitian matrix) auxiliary field a which imposes a unitarity constraint on U{x) 



Z = J DUDaDAi exp{~T[U,A,a]) (4.25) 



M 



r[f/. A, a] = ^^^tr I d'x [{d, - lAi) U\x) (9, + lA^) U{x) + a {u^U - l) 



^The energy, of course never becomes negative, since equation (4.22) contains a divergent M 



independent piece. Here we consentrate only on the M-dependence of E. 
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The role of the vector field Ai is to impose a U{1) gauge invariance, and thereby to 
eliminate one degree of freedom. As far as the thermodynamical properties are concerned, 
its effect is only felt as an 0(1/A^^) correction. At the level of accuracy of the mean field 
approximation, we can safely disregard it, which we do in the following. The mean field 
equations are: 

< f/tf/ >= 1 (4.26) 

< at/ >= (4.27) 
From equation (|4.2?| ) it follows that either < a >= 0, < U >^ (the ordered, bro- 



ken symmetry phase with massless Goldstone bosons), or < a >^ 0, < U >= (the 
disordered, unbroken phase with massive excitations). We are mostly interested in the 
disordered phase, since there the mean field approximation should be reliable. Since the 
symmetry is unbroken, the expectation value of a should be proportional to a unit matrix 



< >= cr'lap (4.28) 



Equation (^1261) then becomes 



2^2(M) fM d^k 1 N^g'^{M) a M 



^ ' \ - — -— = — 1 arctan — = (4.29) 

M Jo (27r)3 A;2 + ct2 \ M a J ^ ' 

The gap equation ( |4.29| ) has solution only for couplings (temperatures) g^{M) larger than 
the critical coupling (temperature) g^, which is determined by the condition that cr = 



The low momentum mode contribution to the ground state energy is 



2 



N^M / - — -— = — -M 

Jo 27r 3 + a2 27r2 



-M-^ - a^M + cr^ arctan — 
3 a 



(4.31) 



The final mean field expression for the ground state energy density is (we do not distin- 
guish between A^^ and A^^ — 1 since we have neglected the contribution of the U{1) gauge 
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field - the errors are of order and are definitely smaller than the error introduced 

by using the mean field approximation in the first place) 



2 ac 



(4.32) 



15 M^a{M)_ 

where a(M) is the QCD coupling at the scale M, ac is given by equation ( [4.30|) , and a 
is determined by 

— arctan — — ) J^c (4.33) 
M a a{M) ^ ' 

The energy as a function of M is plotted on Fig.l for = 3. Qualitatively it is the 
same for any A^. The minimum of the energy is obviously at the point q;(M) = ac- Using 
the one loop Yang - Mills /5 function and Kqcd = 150Mev, we find for = 3 

24 

M = AqcdC" = 8.86Aqcd = 1.33Gev (4.34) 

To see, what is the phenomenological significance of this number we have calculated 
the value of the gluon condensate (a/n) < F^^F^^ >= {2a /n) (< ^ > - < ^ >). 
After some calculations which are straightforward and do not contain any new ingredients, 
we get[| 

< 2 L^jv^M^ < 2 >^ 3— AT^M^ (4.35) 

and finally obtain: 

r < K.^::. >= NM'l (-7^ + 7^) = T^M' = 0.008 Gev' (4.36) 



TT 2 V 407r2 247r2y 1207r 

The best phenomenological value of this condensate, is 0.012 Gev^ ||^. Considering, 

that < F2 > is proportional to the fourth power of M, our result is very reasonable. For 

^We have again kept only the M - dependent pieces. Each one of the quantities < Ef" 2 > 
and < i?f 2 > ig of course positive, due to positive UV divergent, but M independent pieces. It 
is easy to check that the energy density E = l/2(< Ef > + < Bf ^ >) = -(l/307r2)A^2^4 



coincides with the first term in equation (4.32), as it must 
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example, changing M by only 10%, from 1.33 Gev to 1.46 Gev would give {a/ii) < F"^ >= 
0.0116 Gev^, in perfect agreement with [^. 

Note that for = 3, the value of the QCD coupling constant at the variational scale 
is ac = 0.26. It is reasonably small, so that the consistency condition for the perturbative 
integration of the high momentum modes is satisfied. However, it is not so small that 
higher order corrections be negligible. We expect therefore that including higher orders 
in perturbation theory can give corrections to our result for a{M) of order 25%. Since 
M depends exponentially on a{M), such change in a may change the value of M by 
a factor of 2 — 3. Consequently our result for F"^ should be taken only as an order of 
magnitude estimate. This is usually the case in theories with logarithmically running 
coupling constants. The best accuracy is always achieved for dimensionless quantities, 
since those usually are slowly varying functions of a. The overall scale depends on a 
exponentially, and therefore always has the largest error. It would be interesting to 
calculate some dimensionless quantities, such as the ratio of the square of the string 
tension to the SVZ condensate, in our approach p!0[] . 

Another uncertainty comes from the use of the mean field approximation. As a rule, 
mean field approximation gives a good estimate of the critical temperature. Sometimes, 
however it gives wrong predictions for the oder of the phase transition. We believe that 
this is indeed the case here. Our results would indicate that the phase transition is second 
order. The mass gap in the sigma model vanishes continuously at the critical point. The 
universality class describing the symmetry breaking pattern {SU{N) ® SU{N)) /SU{N) 
was considered in the context of finite temperature chiral phase transition in QCD. The 



results of e - expansion |]T2| and also numerical simulations |T3[ strongly suggest that the 
phase transition is of first order. In our case there is an additional Zn symmetry in the 
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game. However, if anything, we believe that its presence should increase the latent heat 
rather than turn the transition into a second order one. The reason is, that the gauge 
invariant theory allows existence of topological defects - the strings, and condensation 
of topological defects frequently leads to discontinuous phase transitions. 

Nevertheless we believe that our results are robust against this uncertainty. The mean 
field approximation should be reliable in the regime where the mass gap in the sigma 
model is not too small. At the point M = i.bAgcD we find 

a = 0.23M, a{M) = 0.38 (4.37) 

Since the gap is of the order of the UV cutoff, the mean field approximation is reliable in 
the vicinity of this point. The perturbation theory is also still reasonable at this value of 
a. The fact, that the energy is negative and has a minimum for some a{M) < .38, seems 
to be therefore unambiguous. 

We now want to argue, that independently of the mean field calculation, it is physically 
very plausible that the energy is minimized precisely at the critical temperature, on the 
disordered side of the phase transition (if it is of the first order). Consider first, the 
contribution of the high momentum modes to the ground state energy, equation ( [4.22|) . 
It is proportional to with a fixed (M - independent) proportionality coefficient x = 
{N'^ — l)/107r^. Consider now the low momentum contribution in the large M region 
equation ( [4.23| ). It is again proportional to with the coefficient ?/o = {^"^ ~ l)127r^. 
The proportionality coefficient of the low momentum contribution at the phase transition 
point, according to our calculation is twice as big yc = 2A^^/127r^ (we disregard the 
difference between N"^ and A^^ — 1). This is very easy to understand physically. In the large 
M - low temperature regime the global symmetry of the sigma model SU (N) ® SU (N) 
is broken down spontaneously to SU{N). This leads to appearance of A^^ — 1 massless 
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Goldstone bosons. In fact, at zero temperature, those are the only propagating degrees of 
freedom in the modeL All the rest have masses of the order of UV cutoff, and therefore do 
not give any contribution to the internal energy. Now, when the temperature is raised (M 
is lowered), the Goldstone bosons remain massless and other excitations become lighter. 
If the transition is second order, at the phase transition point the symmetry is restored, 
and one should have a complete multiplet of the SU {N) (g) SU {N) symmetry of massless 
particles. The dimensionality of this multiplet is 2{N'^ ~ !)• The contribution of every 
degree of freedom to the internal energy is still roughly the same as at zero temperature. 
This is so, since, although at the phase transition the particles are interacting, critical 
exponents of scalar theories in 3 dimensions are generally very close to their values in a 



free theory The internal energy at this point therefore should be roughly twice its 
value at zero temperature. Moving now to higher temperatures, all the particles become 
heavier, and therefore their contribution to internal energy decreases. The internal energy 
therefore should have a maximum at the phase transition temperature. 

Note, that the ground state energy of the Yang - Mills theory, is the difference between 
the high, momentum contributions and the internal energy of the low mode sigma model. 
Already at zero temperature, these two contributions differ only by 20%, and that is why 



the coefficient in the expression equation (|4.24|), even though positive, is so small. At the 



critical point, where the low momentum mode internal energy is twice as large, the chances 
of the slope becoming negative are very good. This is indeed precisely what happens in 
our mean field analysis, but according to the previous argument this in large measure is 
independent of the approximation. If the phase transition is first order one should be more 
careful. The internal energy then changes discontinuously across the phase transition. The 
particles in the disordered phase are always massive, and the internal energy is smaller 
that in the case of the second order phase transition. However, if the transition is only 

27 



weakly first order the same argument still goes through (the fact that the mean field 
predicts second order phase transition, may be an indication that if it is in fact first 
order, it is only weakly so). In fact, it does seem very likely that the ground state energy 
will become negative, since all is need for that, is that the sigma model internal energy 
grows by 20% at the phase transition relative to the zero temperature limit. Moreover, 
in this case there will be a finite latent heat, which means that the internal energy in the 
disordered (high temperature) phase is higher. The ground state energy, therefore, will 
have its minimum in the disordered phase. 

We believe, therefore, that our results are qualitatively correct, and will survive the 
improvement of the approximation. 



5. Wilson loop and area law 

The next interesting question is, whether the variational state we found describes the 
physics of confinement. The relevant quantity to calculate is the Wilson loop 



W{C) =< tr Pexp £ dxiA'^r"^ > (5.1) 



When averaging over A we must take into account the P-ordering of the exponent - the 
simplest way to do it is to introduce new degrees of freedom living on the contour C 



which, after quantization, become the SU{N) matrices r° |15]. We shall consider here 



how it works in the case of SU{2) group - the generalization of this construction to an 



arbitrary Lie group has been discussed in [|T5 



The construction is based on the observation made in [TO that instead of considering 
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the ordered product of matrices one can consider the correlation function 



Dn{t)n''{ti)rf{t2) . . . rf{tk) exp 



i{S + 1/2) / dXe^,^.e'''"'n''^^n''^^n 



where 5* is the spin of representaion, i.e. for fundamental representation S = 1/2, n'^{t) 
is a unit vector ra^ra" = 1 living on a contour C (t is a coordinate on the contour) and S 
is an arbitrary two-dimensional surface with the boundary C = 611. The two-dimensional 
action (here and later we shall concentrate only on case 5* = 1/2) 

S[n] = l^ci^^e^.e'^'^'^aya.n^ (5.3) 

depends only on values n"'[t) at the boundary. The variation of the action is 

55 = £ dte''^''rfdtn^5rf (5.4) 

Here we have used the fact that 5n""n"' = (because n'^n'^ = 1) and thus e^ue"''^'^d^n°'di^n^5n'^ - 
0. It is easy to see that < n"-(ti)n''{t2) . . . n^itk) > depends only on the ordering of ti, . . . 
- as it should. To see this and the fact that n°'{t) behaves effectively as r'^ let us make 
local field reparametrization 

n"(t) ^ n''{t) + e'^^"^]^(t)n"(t) (5.5) 

under which the action variation (|5.4| ) is 5S = — §(jdth°'{t)VL°'{t) and one gets the Ward 
identities (it is important to remember here, that correlators in any QFT are averages of 
the T-ordered products) 

^ < n"(t)n^(ti) . . . n\tk) >= 'Jl^it ' t^^^ < n^{t)n\t,) . . . . . . n'{tk) > (5.6) 

where n'^{ti) means the exclusion of this term from the products of the fields in a correlator. 



From ( |5.6| ) one can conclude immediately that correlation function indeed depends only 
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on ordering of ti, . . . and the following equal time commutaion relations hold 



(5.7) 



which means that one can substitute by a Pauli matrix t°'/2. As a result one can 
represent the Wilson loop (|5.1|) in the form 



W{C) =< J Dn{t)exp ^ j^d^i 



exp (i^ j>^dxiA'^{x{t))n''{t) ) > (5.8) 



and now we can average over Ai using (|3.4|) and 



< exp (ig j!' dxiA°;{x{t))n'^{t)^ >a= exp (y^-W j dxia1{x{t))rf-{t) 

dhdt2X,{h)Ut2)n\h)n\t2){M-^Y\x, y) ] (5.9) 



where a" was defined in ( [4.3|) . Now the Wilson loop can be calculated as the average over 
two scalar fields: U{x) living in the whole space and n"(^) living on a two-dimensional 
surface S such that C = 5S 



W{C) = j DU j Dnexp(-r[f/] +i5[n])exp (^-ig j>^dxia1{x{t))n\t) 

dhdt2X,{ti)Uh)n\ti)n\t2){M-^Y\x,y)] (5.10) 



exp 

V 2JcJc 

In the infrared limit, which is of main interest to us here, we can use (|4.10| ) to simplify 
( |5.10|) and get 

W{C) = J DUlJ Dnexp{-TL[U]+tS[n])exp £ c/x^A^ i(x(t))n"(t)^ 

exp (-^ £ £ dhdt2X,{ti)yi{t2K{h)n''{t2)G{x - y)^ (5.11) 
J DUHexp{-TH[U])exp £dxAl MtYSlVit)^ 

Using equation ([4.20|) one can see that the last term in in ( ^.12|) after integrating over the 
Uh becomes equal to the second term and one gets finally 

WiC) = I Dnexp {tS[n])exp j^dhdt2X,{h)Uh)n\h)n\t2)G{x - y) 
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(5.12) 



where integrating DU is over low-energy modes only and T\U] is the corresponding low- 
energy action. Since G{x — y) is short range, the term 



gives only perimeter dependence and one can neglect it when calculating the string tension. 
Then it can be shown, rewriting n°'{t) as the and performing some simple algebra, that 
the calculation of the Wilson loop is closely related to the calculation of the vacuum 
expectation value of the monodromy operator 



of the sigma model is = SU{N)/Zn, and = Zj^, this factor can take on values 

exp 127171/ N. It has a natural interpretation in terms of the topological defects in the sigma 
model. As mentioned already, the topology allows existence of strings. The string 
creation operator and the operator M satisfy the commutation relations of the t'Hooft 
algebra [|l^. Therefore, in the presence of a string, the operator M has expectation value 
expi27m/N, where n is the linking number between the loop C and the string. As we 
have argued, the sigma model is in the disordered phase. Usually, this means that the 
topological defects are condensed. The vacuum of the sigma model must have therefore a 
large number of strings, and the VEV of Al, probably, will average to zero very quickly, 
and for large loops will have an area law W{C) ~ exp {—a' A). Stricktly speaking, for this 
to happen, one needs not only a large number of strings, but also a large fluctuations in 

this number, but those, usually come together. 

*In fact, the Wilson loop does not reduce to M, but rather to trPexp §q dliU'^ dilJ^ . We 
believe, however, that qualitatively its behaviour should be similar. 




(5.13) 




(5.14) 



in the low momentum sigma model with an effective action r[f/][| Since the target space 
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We also would like to mention that the model of two fields - U and n"- - defined in ( |5.12|) 
is of some interest in itself. For example one can study how nonperturbative fluctuations 
of both flelds - strings and Skyrmions for U and instantons for n interact with each 



other - these questions as well as a calculation of a' will be considered in [|2 • 

An interesting point is, that if one couples fundamental fermions to the Yang Mills 
fields, the effective sigma model will not have a Zn gauge symmetry any more. The origin 
of this Z]\f symmetry is the fact, that the Yang Mills fields do not transform under the 
center of the gauge group. Fundamental fermions, however, do transform nontrivially, 
and therefore the sigma model action will depend on these matrices U. The target space 
now therefore is SU{N), rather than SU{N)/Zn, and is simply connected. The topology 
of the target space does not allow strings any more. Therefore, if it is true, that it is the 
condensation of these objects, that is responsible for the area law for the Wilson loop, the 
area law will disappear. This is in complete agreement with one's expectations, that in a 
theory with fundamental charges, an external test charge can be screened, and therefore 
there is no area law for the Wilson loop. 



6. Discussion and Conclusion 

In this paper we have presented a simple variational calculation of the Yang Mills 
ground state WF. Our trial states preserved gauge invariance explicitly. The results are 
encouraging. We find that the energy is minimal for a state which is different from 
the perturbative vacuum, even though the perturbative vacuum state was included in 
our variational ansatz. Dynamical scale generation takes place and the gluon (SVZ) 
condensate in the best variational state is nonzero. 
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It is interesting to note, that from the point of view of the effective sigma model, the 
energy is minimized in the disordered (unbroken) phase. In other words, the fluctuations 
of the field U are big, unlike in the perturbative regime (high momentum modes), where 
U is very close to a unit matrix. Prom the point of view of the original WF this means 
that the off diagonal contributions, coming from the Gaussian gauge rotated by a slowly 
varying gauge transformations, are large. This is telling us, that it was indeed necessary 
to project the initial Gaussian onto gauge invariant state. Without doing this, we would 
miss the contributions of the off diagonal elements to the energy expectation value. 

There is still a lot of work to be done, even in the framework of our variational ansatz. 
Our present paper should be considered only as an exploratory research. Of course, 
coupling the fermions is a very interesting question in itself. It seems to us, that it should 
be possible to treat a theory with fermions in basically the same variational approach as 
presented here. It would be then very interesting to see the chiral symmetry breaking and 
calculate fermionic condensates. 

Quite apart from this, there are several technical points that can be improved. First, 
we are planning to extend the RG calculation to take into account the one loop contri- 
bution of the high momentum modes. This might require to change a variational ansatz 
slightly. One may have to consider not the gauge projected Gaussians, but gauge pro- 
jected products of Gaussians and polynomials of the fourth order in the fields. This does 
not change the level of complexity of the calculation. 

It would also be desirable to have better methods to deal with the low momentum 
sigma model, especially since we suspect that the mean field approximation does not give 
the correct order of the phase transition. Although we do not expect the variational 
parameter to be very sensitive to this, the vacuum condensates can depend strongly on 
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the mass gap of the sigma model. 

Finally, there is one more direction, in which the calculation can be extended. In this 
paper we have adopted the simplest ansatz for the width of the Gaussian G, based on the 
argument, that it should be short ranged. The Fourier transform of our propagator goes 
to a constant at zero momentum. This, however is no the only possible form of a short 
range correlator. It could have a different small momentum behaviour. It is quite possible 
that the small momentum behaviour is very important. One could therefore introduce 
an additional variational parameter 7, assuming the asymptotic small momentum depen- 
dence of the function G to be of the form . This will only affect the last step of our 
calculation. The action of the effective low momentum model will have extra derivatives. 

In conclusion, it seems to us that the type of the variational approximation presented 
here is manageable, and also gives some preliminary interesting results. It therefore 
warrants further work along the lines described in this section. 
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Figure Caption. 

Fig.l The energy density of a variational state as a function of the variational parame- 
ter M in units of Aq^d- The energy is only shown for M < 8.86Aqc7d, which corresponds 
to the disordered phase of the effective low momentum a model. Close to the phase 
transition point in the ordered phase, the mean field approximation is not applicable. 
Far from the phase transition point, at large M the energy density is a monotonically 
increasing function of M given in equation ( [4.24 ). 
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